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8 . ABSTRACT 

Context. Classical novae (CNe) have been found to represent the major class of supersoft X-ray sources (SSSs) in our neighbour 
galaxy M 3 1 . 

Aims. We determine the properties and evolution of the two first SSSs ever discovered in the M 31 globular cluster (GC) system. 
Methods. We have used XMM-Newton, Chandra and Swift observations of the centre region of M 31 to discover both SSSs and to 
determine their X-ray light curves and spectra. We performed detailed analysis of XMM-Newton EPIC PN spectra of the source in 
Bol 111 (SSI) using blackbody and NLTE white dwarf (WD) atmosphere models. For the SSS in Bol 194 (SS2) we used optical 
monitoring data to search for an optical counterpart. 

Results. Both GC X-ray sources were classified as SSS. We identify SSI with the CN M31N 2007-06b recently discovered in the 
M 31 GC Bol 111. For SS2 we did not find evidence for a recent nova outburst and can only provide useful constraints on the time of 
the outburst of a hypothetical nova. 

Conclusions. The only known CN in a M 31 GC can be identified with the first SSS found in a M 31 GC. We discuss the impact of 
our observations on the nova rate for the M 3 1 GC system. 

Key words. Galaxies: individual: M 31 - novae, cataclysmic variables - stars: individual: Nova M31N 2007-06b - globular clusters, 
individual: Bol 111, Bol 194 - X-rays: galaxies 

x. 

^ , 1. Introduction and CAL 87 ( see IParmar et all 1 1 9971. and references therein), 

both of which are located in the Large Magellanic Cloud. These 
Supersoft X-ray sources (SSSs) are a class of X-ray sources objec ts show eclipses (CAL 87) or rar e X-ray off states (CAL 
that were first characterised based on ROSAT observations (see 83) dKahabka & van den Heuvell Il997h. but are rat her perma- 



e.g. iGreiner et al.| [19911). These sources show extremely soft nent SSSs. On the other hand. IPietsch et all (T2005^h found out 

X-ray spectra, with little or noradiation at energies above 1 that classica i novae (CNe) represent the major class of SSSs 

keV (see e.g. [ Parmar et alj |j_998j), that can be described by in ou r neighbour galaxy M 31 (distance 780 kpc, iHollandl 

blackbody temperatures typically in the range of 15-80 eV [19981: IStanek & Garnavichl[l998l CNe can appear as luminous 

(see | Kahabka & van den Heuvel| | 1997i and references therein). (L ^ 1Q 38 erg g -ij trans i e nt SSSs that seem to go thr ough a single 

SSSs were originally believed to be hydrogen burning white outburst that can last from months to several years (IPietsch et al.1 

dwarfs (WDs) in binary systems, where the WD steadily |2QQ7bb 
burns hydrogen rich matter accre ted from its companion star 

(Kahabka & van den Heuvell9 97). However, the class of SSSs Classical novae (CNe) are thermonuclear explosions on the 
is quite inhomogen eous. Prototypical sources are, on the one surface of WDs in cataclysmic binaries that result from the trans- 
hand, CAL 83 (see IParmar et"al]|l998L and references therein) fer of matter from the companion star to the WD. The trans- 
ferred hydrogen-rich matter accumulates on the surface of the 
* Partly based on observations with XMM-Newton, an ESA Science WD until hydrogen ignition starts a thermonuclear runaway in 
Mission with instruments and contributions directly funded by ESA the degenerate matter of the WD envelope. The resulting expan- 
Member States and NASA sion of the hot envelope can cause the brightness of the WD 
** current address: California Institute of Technology, Pasadena, CA to rise by more than nine magnitudes within a few days, and 
91 125, USA mass to be ejected at high velocities (see lHernanjl2005l I Warner! 
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I1995L and references therein). However, a fraction of the hot en- 
velope can remain in steady hydrogen burning on the surface 
of the WD (IStarrfield et al.l 1 974. [Sala & Hernanz.20051) . power- 
ing a supersoft X-ray source that can be observed dir ectly once 
the ejected envelope becomes sufficiently transparent (IStarrfieldl 
11989. lKrauttet1l2002l) . 

The duration of the supersoft phase is related to the amount 
of H-rich matter that is not ejected and also depends on the lu- 
minosity of the white dwarf. More massive WDs need to accrete 
less matter to initiate the t hermonuclear runaway , because of 
their higher surface gravity dJose & Hernanzl[l998l) . In general, 
more massive WDs retain less accreted mass after the explosion, 
although this also depends on the accretion rate, and reach larger 
luminosities (lYaron et al.ll2005l) . Thus, the dur ation of the SSS 
state is inversely related to the m ass of the WD dSala & Hernanzj 
120051: iTuchman & TrurarJll998l) . In turn, the transparency re- 
quirement mentioned above implies that the time of appear- 
ance of the SS S is determined by the fraction of mass ejected 
in the outburst (Ha chisu & Ka to 2006). X-ray lightcurves there- 
fore provide important clues on the physics of the nova outburst, 
addressing the key question whether a WD accumulates matter 
over time to become a potential progenitor for a type la super- 
nova (SN-Ia). 

Due to its proximity to the Galaxy and its moderate Galacti c 
foreground absorption (N n = 0.7xl0 21 cm~ 2 . IStark et al.lll992T) 
M 31 is a unique target for CN surveys. Starting with iHubbld 
(1929|) the majori ty of these surveys that were conducted in 
the past (see e.g. iHenze et all 120 08b, and references therein) 
classified optical transients as novae just by their lightcurve. 
Following the steep rise, the luminosity of different CNe declines 
with different speed, which allows a pheno menological classifi- 
cation of these objects by their speed class dPavne-Gapos chkin 
Il964l) . Eventually, novae in M 31 will fade back to invisibility, 
since the binary systems in quiescence are too faint to be ob- 
served at the distance of M 31. Only recently, nova monitoring 
programs for M 3 1 were established, that include fast data anal- 
ysis and that therefore provi de the possibility to c onduct follow- 
up sp ectroscopy (see e.g. IHenze et all l2008at iDi Mille et al.1 
2008) and to classify CNe within the system of lWilliamsl (I1992TT 

Almost all optical surveys for CNe in M 3 1 that were con- 
ducted in the past searched for suddenly appearing objects that 
have not been visible before and fade back to invisibility in days 
to weeks. This condition is certainly not fulfilled by CNe in rela- 
tively bright GCs, where the optical background light of the GC 
itself makes a photometric discovery of a nova outburst much 
more complicated. Therefore, the connection of CNe to SSSs 
in X-rays provides a useful possibility to detect CNe in GCs. 
Novae in GCs are rare. There were just two sources known so 
far that likely fit this definition. One of them was seen in the 
Galactic GC M 80 wh ereas the second n ova was found in a GC of 
the galaxy M 87 (seelShara et al. 2004, and references therein). 
According to lShara et al.l (120041) a third candidate (nova 1938 in 
the galactic GC M 14) is less likely to be a GC nova. Recently, 
Shafter & Ouimby (2007) reported the very first nova found in a 
M 3 1 GC (M3 IN 2007-06b). 

Independent from their connection to CNe, SSSs in GCs 
are rare objects. There was just one SSS known in GC up 
to now: the transie nt IE 1339 8+283 7 in the Galactic GC 
M3 (NGC 5272) dDotani et al.l [1999. IVerbunt et al.l fT995l) . 
IE 1339.8+2837 was first discovered with the Einstein satel- 
lite with a hard X-ray spectrum and a low luminosity of 
4 xlO 33 ergs" 1 dHertz & Grindlavl [19831) but showed a very 
soft high lu minosity state dur i ng 1991-92 when observed 
with ROSAT dDotani et al.lll999l) . iDotani et all d 19991) give for 



IE 1339.8+2837 a blackbody temperature of kT ~ 36 eV and a 
luminosity of ~ 10 35 erg s _1 , which is significantly lower than 
the observed peak luminosities of other SSS, including the two 
sources discussed in this work. They discuss this object as a cat- 
aclysmic variable (CV) system that may be a dwarf nova in- 
cluding a mass i ve WD . The CV interpretation is supported by 
Edm onds et al.l d2004l) who used Hubble Space Telescope obser- 
vations to discover an optical counterpart to IE 1339.8+2837. 
They suggest that magnetically channeled accretion could ex- 
plain the peculiarities of IE 1339.8+2837. 

Recently, we rep orted, based on preliminary data analysis 
dPietsch et al.ll2007al) . the very first SSS (hereafter SSI) that was 
found in a M 31 GC. During the same observing run where we 
found this source we detected a second luminous a nd previously 
unkno wn SSS (hereafter SS2) in another M 31 GC dHaberl et all 
l2007allbl) . Both sources are transients and will be examined in 
detail in this work. 

The structure of the paper is as follows: In Sect. 2 we de- 
scribe our X-ray observations and the properties of the two SSSs. 
The optical data available for Bol 194 and the data analysis are 
presented and interpreted in Sect. 3. Finally, in Sect. 4 we sum- 
marise our results and discuss the impact that two GC novae 
within one year would have on the nova rate for the M 3 1 GC 
system. 



2. X-ray observations and data analysis 

In the context of the XMM-Newton/C/ianiira M 3 1 nova mon- 
itoring projecfl we obtained five 20 ks Chandra HRC-I ob- 
servations of M31 starting in November 2007. We detected 
two new sources in the M31 GCs Bol 111 and Bol 194 
dGalleti et al.ll2004l) . For both sources follow-up ToO observa- 
tions with the Swift X-r a y Telescope (XRT) were requested 
(e.g. Piet sch et ID l2007at (Kong & Pi Stefanol [2007b) as well 
as additional Chand ra ACIS-S observations by another group 
dGalache et al.ll2007l) . We also used two observations obtained 
within the XMM-Newton M 31 large survey project- Details 
on the observations used are given in TablesQ] and [2] for the 
sources in Bol 111 and Bol 194, respectively. The tables list 
the telescopes and instruments used, the observation identifica- 
tions (ObsIDs), the dates as well as source count rates and lumi- 
nosities or upper limits, respectively. X-ray light curves of both 
sources are shown in Fig.Q] Figure[2] shows the fields of view 
of Chandra and XMM-Newton observations with representative 
positions. 

We analysed all available observations using mission de- 
pendent source detection software as well as the HEAsoft 
package v6.3, including the spectral analysis software XSPEC 
vl2.3.1. In all our XSPEC models we used the Tubingen- 
Boulder ISM absorption (TBabs in XSPEC) model to- 
gether with the photoelectric absorption cross-sections from 
iBalucinska-Church & McCammonl d 19921) and ISM abundances 
from [Wilms et al.l (12000) . For the individual telescopes we ap- 
plied the following data reduction techniques: 

The Swift XRT data were analysed using the HEAsoft 
XIMAGE package (version 4.4) with the sosta command 
(source statistics) for estimations of count rates. We took into 
account the XRT PSF of the sources that we computed with 



1 http://www.mpe.mpg.de/~m3 lnovae/xray/index.php 

2 XMM-Newton Large Program: The X-ray source population of the 
Andromeda galaxy M 31; PI: W. Pietsch 

see http://www.mpe.mpg.de/xray/research/normal_galaxies/m31/lp.php 
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the command psf, as well as exposure maps that were cre- 
ated with the XRT software task xrtexpomap within XIMAGE 
Astrometry was done using the HEAsoft routine xrtcentroid 

For the XMM-Newton data we applied a background screen- 
ing and used the XMMSAS v6.6 tasks eboxdetect anc 
emldetect to detect sources in the image and perform astrome- 
try and photometry. For computing upper limits we added an ar- 
tificial detection at the position of the source to the eboxdetect 
list. This list was used as input for an emldetect run (with 
fixed positions and likelihood threshold of zero) that derivec 
the observed flux and upper limit for all objects in the list 
We give 3 times the background flux as the 3cr upper limit 
To obtain astrometric ally-corrected positions we selected opti- 
cal sources, with prope r-motion corrected positions, from the 
USNO-B1.0 catalogue dMonet et alj 120031) . and checked thai 
only one optical source was in the error circle of the cor- 
responding X-ray source. We only accepted sources correlat- 
ing with globular clus ters from the Revised Bo l ogna Catal ogue 
(V.3.4, January 2008; iGalleti etal.ll2004l 120051 120061 120071) or 
with fore ground stars, characteri zed by their optical to X-ray 
flux ratio dMaccacaro et al.l H988) and their hardness ratio (see 
IStiele et al.ll2008l) . We then used the XMMSAS task eposcorr 
to derive the offset of the X-ray aspect solution. 

We reduced the Chandra observations with the CIAO v3.4 
(Chandra Interactive Analysis of Observations) software pack- 
age. The source detection was done with the CIAO tool 
wave detect. Similar to the procedure described by Elvis et all 
(2009) for the Chandra COSMOS survey, an adapted version of 
the XMMSAS tool emldetect was used to estimate background 
and exposure corrected fluxes and count rates for the detected 
sources. Both sources are located near the edge of the HRC- 
I field of view and therefore both, photometry and astrometry, 
suffer from relatively large errors. 
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Fig.l. X-ray light curves of SSI and SS2 (0.2 - 1.0 keV) 
obtained from XMM-Newton (black), Chandra HRC-I (red), 
Chandra ACIS-S (blue), and Swift XRT (green) (see also 
Tables [Tl2"l ). The upper limit for SSI is indicated by a down- 
pointing arrow, the upper limits for SS2 are marked by open 
lozenges. An open square indicates the lower limit luminosity 
for SSI on Chandra ObsID 8526. Crosses symbolise luminosi- 
ties of detections. Error bars are not shown, since the (statistical) 
errors are mostly smaller than the size o f the symbols. There are 
no errors given bv lGalache et alj (|2007) for the ACIS-S observa- 
tions. The vertical dotted line indicates the d ay of the first detec- 
tion of nova M31N 2007-06b in the optical dShafter & Ouimbvl 
120071) . 



2.1. Supersoft source in Bol 111 

In our first Chandra observation of the A06 monitoring cam- 
paign, starting on 2007-11-07.64 UT (ObsID 8526), we de- 
tected a new source (SSI) at the very edge of the HRC-I field 
of view (Piet sch et al.l|20 07a). SSI remained active during the 
following Chandra monitoring observations (ObsIDs 8527-30) 
in 2007 November and December. We followed the light curve 
of SSI with Swift ToO observations (ObsIDs 00031017001-12, 
2007 November - 2008 January) and also found the source to 
be still visible in our XMM-Newton M 31 monitoring obser- 
vations (ObsIDs 0511380201 and 0511380601, 2008 January 
- February). In TableQ] we present details on all observations. 
Due to the location of SSI near the edge of the Chandra HRC- 
I field of view during all observations (see Fig. [2] for the loca- 
tion of the Chandra field and the source), we used the XMM- 
Newton observations to perform precise astrometry. The offset- 
corrected XMM-Newton position of the source was determined 
to be RA(J2000) 00:42:33.21, Dec(J2000) +41:00:26.1 with a 
3<t error of 1 '.'6, including the uncertainty of the offset correc- 
tion. These coordinates are in good agreement (distance = 0'.'5), 
within the errors, w ith the position of the M31 GC Bol 111 
(IGalleti et al l 12001 00:42:33.16, +41:00:26.1). Therefore, we 
assume that SS 1 is situated within the globular cluster. Note, that 
no X-ray source was previously known in this GC. See Fig.[3]for 
a visualisation of the agreement of optical (DSS POSS-II Red) 
and X-ray positions. 

To perform spectral analysis of SSI we used XMM- 
Newton observations obtained on 2008-01-05.99 UT and 2008- 
02-09.31 UT (ObsIDs 0511380201 and 0511380601). We ex- 
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Fig. 2. ROTSE-III M 31 central field (see Sec.©. Overlaid are 
the four Super-LOTIS fields (black squares), the Chandra HRC- 
I field for ObsID 8527 (big blue square) and the XMM- 
Newton field for ObsID 051 1380201 (big red square). Indicated 
are the position of Bol 111 (red circle) and Bol 194 (green 
square). 
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Table 1. X-ray observations of SSI in Bol 1 1 1. 



Telescope/Instrument" 


ObsID 


Exp. time* 


Date r 


Offset'' 




Count Rate c 




Lo.2-1.0 ' 






[ks] 


[UT] 


[d] 




[ct s '] 




[erg s -1 ] 


XMM-Newton EPIC PN 


0505760201 


49.2 


2007-07-22.55 


33 




< 1.1 xlO 3 




<1.4 xlO 36 


Chandra HRC-I 


8526 


18.7 


2007-11-07.64 


141 


> (19 


± 0.4) xlO 2 


>(5.3 


± 1.1) xlO 38 


Chandra HRC-I 


8527 


20.0 


2007-11-17.76 


151 


(2 5 


± 0.2) xlO" 2 


(7.0 


±0.7) xlO 38 


Chandra HRC-I 


8528 


20.0 


2007-11-28.79 


162 


(2.2 


± 0.2) xlO" 2 


(6.0 


±0.6) xlO 38 


Chandra HRC-I 


8529 


18.9 


2007-12-07.57 


171 


(2.8 


± 0.3) xlO" 2 


(7.7 


± 0.7) xlO 38 


Chandra HRC-I 


8530 


19.9 


2007-12-17.49 


181 


(2 5 


+ 3} xl0~ 2 


(7.1 


± 0.7) xlO 38 


Swift XRT 


00031017001/2 


7.1 


2007-11-18.40 


152 


(1.18 ± 


- 0.15) xl0~ 2 


(11.6 ± 1.4) xlO 38 


Swift XRT 


00031017003 


3.0 


2007-11-13.02 


177 


(0.8 


± 0.2) xlO 2 


(8.1 


± 1.9) xlO 38 


Swift AK1 


0003101 /004 


3.0 


200 /-1 2- 14.02 


178 


(1.1 


i A 1\ w1 A— 2 

± 0.2) XlO 


(10.6 


± 2.2) xlO 38 


Swift XRT 


00031017005 


3.2 


2007-12-15.03 


179 


d 1 


± 0.2) xlO" 2 


(11.1 


± 2.1) xlO 38 


Swift XRT 


00031017006 


2.2 


2007-12-20 25 


184 


CI 1 


+ 3} x 1 0~ 2 


(11.1 


± 2.6) xlO 38 


Swift XRT 


00031017007 


2.1 


2007-12-22.39 


186 


(0.5 


± 0.2) xl0~ 2 


(5.1 


± 1.9) xlO 38 


Swift XRT 


00031017008 


2.3 


2007-12-24.33 


188 


(0.9 


±0.2) xlO" 2 


(9.1 


± 2.4) xlO 38 


Swift XRT 


00031017009 


2.3 


2007-12-30.15 


194 


(1.0 


±0.2) xlO" 2 


(10.0 


± 2.4) xlO 38 


Swift XRT 


00031017010 


2.0 


2008-01-03.44 


198 


(1.2 


±0.3) xlO" 2 


(12.2 


± 2.8) xlO 38 


Swift XRT 


00031017011 


1.9 


2008-01-06.25 


201 


(0.4 


±0.3) xlO" 2 


(4.0 


± 2.8) xlO 38 


Swift XRT 


00031017012 


1.7 


2008-01-10.00 


205 


(1.1 


± 0.3) xlO" 2 


(11.1 


± 1.3) xlO 38 


XMM-Newton EPIC PN 


0511380201 


23.0 


2008-01-05.99 


200 


(8.2 


± 0.2) xl0~ 2 


(10.6 


±0.2) xlO 38 


XMM-Newton EPIC PN 


0511380601 


24.0 


2008-02-09.31 


235 


(7.5 


±0.2) xlO" 2 


(12.2 


±0.4) xlO 38 


Swift XRT 


00037718001 


4.8 


2008-05-26.29 


342 


(0.5 


±0.1) xlO" 2 


(4.6 


± 1.3) xlO 38 


XMM-Newton EPIC PN 


0560180101 


17.4 


2008-07-18.26 


395 


(3.0 


± 0.2) xlO" 2 


(8.9 


±0.5) xlO 38 



Notes: ": Telescope and instrument used for observation. 

b : Dead time corrected exposure time of the observation. 

c : Start date of the observation. 

d : Time in days after the discovery of nova M31N 2007-06b in the optical dShafter & Ouimby||2007h on 2007 June 19.38. 
(JD = 2454271). 

e : Source count rates, X-ray luminosities (unabsorbed, blackbody fit, 0.2 - 1.0 keV) and upper limits were estimated according 
to Sect. |2.1| For Chandra ObsID 8526 the source is right on the detector edge, therefore we give lower luminosity limits. 




(a) Bol 111 (b) Bol 194 



Fig. 3. X-ray position error circles for SSI and SS2 as de- 
scribed in Sect. 12. II and Sect. 12. 21 with respect to Bol 111 and 
Bol 194. Blue/big: Chandra HRC-I, red/medium: Swift XRT, 
green/small: XMM-Newton PN. Underlying optical image: 5' x 
5' DSS POSS-II Red. 



tracted the spectra of SS 1 from both observations and fitted them 
simultaneously in order to increase statistics. We used data from 
XMM-Newton's PN detector because of the better sensitivity of 
PN in the soft band compared to both MOS detectors. The spec- 
tra were extracted using XMMSAS task evselect. For both 
spectra only single-pixel events (PATTERN = 0) were selected. 

We fitted the spectra in XSPEC using an absorbed black- 
body approach. Temperature and foreground A^h were both as- 



sumed to be the same during the two observations and only the 
respective normalisations were allowed to vary independently 
from each other. The relative stability of the spectral parame- 
ters is confirmed by a XMM-Newton observation on 2008-07- 
18.26 UT (160 days later than 051 138601), from which we ex- 
tracted a spectrum of SSI that can be fitted by a model with 
similar parameters (see below). The blackbody approach yields 
an acceptable = 1-39 for the best fit values of kT = 48^ 
eV and N H = 2.3 + 0.1 xlO 21 cirr 2 . This fit is shown in Fig.g] 
and the associated contour plot is given in Fig. [5] We computed 
the unabsorbed EPIC PN X-ray luminosities, in the range 0.2 - 
1.0 keV, from the best fit model in XSPEC and used the best 
fit values to create in XSPEC fake spectra (command fake it) 
to infer the energy conversion factors (ecf) for the Swift XRT 
(ecfyfir) and the Chandra HRC-I (ecf hrc-i)- The ecf values are 
given in Table|3] and were used to convert our Swift XRT and 
Chandra HRC-I count rates to unabsorbed lumi nosities. All lu- 
minosities assume a distan ce to M 31 of 780 kpc (Holland 1998; 
Stanek & Gar navichlll998l) and are presented in TableQ] Note, 
that for our first detection of SSI in the Chandra HRC-I obser- 
vation 8526 no photometry is possible due to the location of the 
source on the edge of the detector. The blackbody fit parameters 
and derived values, like luminosities for ObsID 051 1380201, are 
given in Table[3] 

The unabsorbed X-ray luminosities inferred from the black- 
body fit are in the order of 10 39 erg s 1 , and therefore signif- 
icantly exceed the Eddington luminosity of the hydrogen rich 
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Telescope/Instrument 


ObsID 


Exp. time^ 


Date 1 ' 


Offset'' 




Count Rate e 




Luminosity^ 








TUTl 

[kj i j 


L U J 




let s" l l 






XMM-Newton EPIC PN 


0505760201 


49.2 


2007-07-22.55 







< 2.9 xlO" 3 




< 1.6 xlO 36 


L-nanara /M^lo-o 




J.U 


inm i i ni 1 C 
ZUU / - 1 1-Uj. 1 o 


i n A 

1 U4 




LA X 1U 




j. 4 X1U 


Chandra HRC-I 


8526 


18.7 


2007-11-07.64 


108 


(3.55 ± 


0.25) xlO" 2 


(8.41 


±0.59) xlO 37 


Chandra HRC-I 


8527 


20.0 


2007-11-18.76 


118 


(2.95 ± 


0.25) xlO" 2 


(7.00 


±0.60) xlO 37 


Swift XRT 


00031027001 


7.3 


2007-11-24.34 


125 


(8.0 


± 1.2) xlO" 3 


(6.35 


± 1.04) xlO 37 


Chandra ACIS-S 


8187 


5.0 


2007-11-27.16 


128 




< 6 xlO" 4 




< 9.7 xlO 35 


Swift XRT 


00031027002/3 


4.7 


2007-12-02.64 


133 


(3.2 


t 1.0) xlO" 3 


(3.92 


± 1.22) xlO 37 


Swift XRT 


00031027004 


3.9 


2007-12-16.77 


147 




< 3.2 xlO- 3 




< 3.9 xlO 37 


Swift XRT 


00031027005 


4.0 


2007-12-30.02 


160 




< 2.0 xlO" 3 




< 2.4 xlO 37 


XMM-Newton EPIC PN 


0511380201 


22.8 


2008-01-05.99 


167 




< 5.1 xlO" 3 




< 3.5 xlO 36 



Notes: Telescope and instrument used for observation. ACIS-S count rates are from lGalache et all d2007h . 
: Dead time corrected exposure time of the observation. 
c : Start date of the observation. 

d : Time in days after the last X-ray non-detection of SS2 on 2007 July 22.55 (JD = 2454304). 
e : Source count rates, X-ray luminosities (unabsorbed, blackbody fit, 0.2 - 1.0 keV) and upper limits were estimated 
according to Sect. 12.21 




Channel Energy (keV) N H (10 21 cm 2 ) 



Fig. 4. XMM-Newton EPIC PN spectra of SSI (crosses) from 
observations 0511380201 (black) and 0511380601 (blue) fitted 
with an absorbed blackbody (solid lines). 



atmosphere of a WD: L Edd = 1.3 X 10 38 ^jerg s" 1 . The 

fact that blackbody fits to SSS spectra produce in general too 
high values of Nh and too low temperatur es, and therefore too 
high luminosities, is well known (see e.g. iGreiner et al.lll99lt 
iKahabka & van den Heuvellll997l and references therein), thus 
the values given in TableQ] define upper limits on the actual lu- 
minosity of SSI. 

Blackbody fits are very simple approximations for SSS emis- 
sion of novae and allow to compare parameters like effective 
temperature for different novae. However, as mentioned above, 
these fits are physically not realistic. Therefore, we tested fitting 
our low -resolution spectra with WD atmosphere models, that are 
based on more physical assumptions. We used a grid of syn- 
thetic ionizing spectra for hot compact s tars from NLTE model 
atmospheres computed by Rauch (2003). These NLTE models 
are plane-parallel, in hydrostatic an d radiative e quilibrium and 
contain all elements from H to Ca dRauch|[l997l) . The models 



Fig. 5. Column density (Nh ) - temperature (kT) contours in- 
ferred from the blackbody fit to the XMM-Newton EPIC PN 
spectra of SSI (see Hig-SJ. Normalisation has been adjusted. 
Indicated are the formal best fit parameters (cross) and the lines 
of constant X-ray luminosity (0.2-1.0 keV, dotted lines). 

were co mputed using the Tubinge n Model- Atmosphere Package 
(TMAP, iRauch & Deetjenl [2003). Elemental abundances are 
fixed to either galactic halo ([X]j = [Y] = 0, [Z] = -1) or so- 
lar ([X] = [Y] = [Z] = 0) ratios. The grids of model atmosphere 
fluxes, as well as FITS tables which can be used in XSPEC, are 
available on-linq^. The XSPEC tables contain temperatures and 
fluxes (binned to 0. 1 A intervals) for fixed surface gravity (logg) 
and elemental abundances. In our case, the available grid param- 
eter space was restricted to models with logg = 9.0, since only 
these tables include temperatures high enough to fit our spectra. 
This restriction clearly limits the significance of our best fits, 
since models with different surface gravity may have provided 
equally good fits with different parameters. Also note, that the 

3 [] : log(abundance/solar abundance) 

4 http://astro.uni-tuebingen.de/~rauch/ 
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assumptions of plane parallel and static are not physically real- 
istic for a nova atmosphere. Furthermore, the spectral analysis 
is limited by the low energy resolution of the EPIC PN spectra, 
which are, due to the faintness of the source, the only available 
spectra with sufficient signal to noise ratio. Great caution should 
therefore be applied when interpreting the results of our test fits. 

We fitted the EPIC PN spectra with both halo and solar 
abundance ratios. Galactic halo abundances are more represen- 
tative of the metallicity of the M 31 GCs dBarmbv & Huchral 
2000). Figure[6] gives the fit for halo abundances and shows that 
the atmosphere model fits the spectrum well only up to an en- 
ergy of about 600 eV (x 1 doL 1 = 1.01) but predicts too high 
flux at energies above, which is clearly indicated by the resid- 
uals. This is true for solar and halo abundances. However, this 
is no surprise since recent detailed X-ray spectroscopy of the 
SSS phase of galactic novae shows that these spectra are com- 
plex (see e.g. iNess et aill2007l who used spectra obtained with 
the XMM-Newton Reflection Grating Spectrometer (RGS) and 
the Chandra Low Energy Transmission Gratings Spectrometer 
(LETGS)). Spectral models currently available may be not suit- 
able for a correct fitting of SSS spectra of novae, especially with 
respect to the elemental abundances of a hot post-nova atmo- 
sphere. 

In Table[3] we compare the parameters of the WD atmo- 
sphere model fits (fitted to the 0.2 - 0.6 keV range) to the pa- 
rameters of the blackbody fit described above. We give unab- 
sorbed X-ray luminosities(L A ), that were computed for the range 
0.2 - 1.0 keV in XSPEC, and bolometric luminosities (L\, \) as 
well as WD radii (R), inferred from the bolometric luminosi- 
ties. The bolometric luminosity for the blackbody fit was directly 
computed from the normalisation of the spectra in XSPEC. For 
the WD atmosphere model we used an absorbed blackbody fit 
with the same temperature and foreground Nh to compute bolo- 
metric luminosity and WD radius. These values overestimate 
the actual luminosity and radius and are therefore upper limits. 
Energy conversion factors (ecf) for the Swift XRT (ecf xrt) and 
the Chandra HRC-I (ecf hrc-i) detectors were again computed in 
XSPEC, using fake spectra as for the blackbody fit, and are listed 
in Table[3] together with the ecf for the XMM-Newton EPIC PN 
(ecf/w). 

The comparison of the fits shows that the WD atmosphere 
models give physically more plausible results. The high value 
of Nh which was required for the blackbody fit (three times the 
foreground absorption) is significantly reduced to a value which 
is just slightly above the foreground level of Nh = 0.7xl0 21 
cm -2 . As a consequence, the super-Eddington luminosity de- 
rived from the blackbody fits is reduced to more realistic val- 
ues for the atmosphere models. The connection of Nh and effec- 
tive temperature can be seen in the confidence contours shown 
in Fig. [5] The effective temperature for the WD atmosphere is 
higher than the blackbody temperature and is revealing the H- 
burning hot layer. Still, one has to be cautious interpreting these 
models, because the elemental abundances are fixed and may 
be not realistic to describe abundances in a burning WD atmo- 
sphere. However, the consideration of model substructures, like 
absorption edges, already causes a reduction of the /Vh necessary 
to adjust the model to the spectra. Truly physical models for su- 
persoft emission from novae are strongly needed and should be 
validated using high resolution spectra obtained with current and 
future missions. 

We checked recent Swift (ObsID 00037718001, starting at 
2008-05-26.29 UT) and XMM-Newton (ToO 0560180101, start- 
ing at 2008-07-1 8.26 UT) data and found that SSI is still visible 
in both observations (see also TableQ] for details). Preliminary 



Table 3. Comparison of SSI spectral best fit parameters and 
derived parameters for blackbody and WD atmosphere models 
with halo and solar abundances. 



1V1UUC1 

(energy range [keV]) 


(0.2 - 0.8) 


WD hnln 
vv lj uaiu 

(0.2 - 0.6) 


WD tinier 
VV LJ hULal 

(0.2 - 0.6) 


K 1 V ) 




61 + 1 

Ul X 1 


7fl + 1 


A'h (lu cm ) 


1.5 ± U.l 


l.U ± U.z 


l.U ± U.z 


xl 


1.39 


1.01 


1.06 


dof 


53 


44 


48 


L, (10 38 erg s"') 


10.6 + 0.2 


1.00 ±0.02 


0.88 ± 0.02 


L bol (10 38 erg s- 1 ) 


2S.T_° j 


2.7 + 0.1 


1.5 ±0.1 


(10* cm) 


7 +L6 

'• u -0.7 


1.21 ±0.07 


0.70 ± 0.04 


ecfpjv (ct cm 2 erg -1 ) 


5.0xl0 9 


5.5xl0 10 


6.1xl0 10 


ecf hrc-i (ct cm 2 erg" 1 ) 


2.6x10* 


3.2xl0 10 


3.7x10'° 


ecW (ct cm 2 erg"') 


7.4x10 s 


8.7x10* 


9.9x10* 



Notes: Luminosities and WD radii refer to the XMM-Newton obser- 
vation 051 1380201. The unabsorbed X-ray luminosity L x is for the 0.2 
-1.0 keV range. Bolometric luminosities and WD radii for the WD at- 
mosphere model are upper limits (see Sec. l2.1l for details). 




Fig. 6. XMM-Newton EPIC PN spectra of SSI from observa- 
tions 0511380201 (black) and 0511380601 (blue) fitted with 
an WD atmosphere model with galactic halo abundances (solid 
lines). 

data analysis was published in Piet sch et alj d2008albl) . For 
Swift we used the blackbody fit described above and in Table[3] 
to convert the XRT count rate to unabsorbed luminosity. For 
XMM-Newton we extracted an EPIC PN spectrum of SSI as 
described above for the two earlier XMM-Newton observations. 
The spectrum can be fitted with an blackbody fit with best fit val- 
ues of kT = 43 ± 8 eV and N H = 2.2+1-2 xlO 21 cnr 2 (y 2 = 1.52). 
These parameters are consistent with the blackbody fit described 
above and result in an unabsorbed X-ray luminosity (range of 0.2 
- 1.0 keV) of (8.9 + 0.5) XlO 38 erg s" 1 . Both observations show 
that the X-ray luminosity of SSI is slowly declining. 

Supersoft spectra like the one of SSI are typical for coun- 
terparts to optical novae and indicate that after the outburst 
there is hydroge n burning going on i n the remaining envelope 
of the WD (see iPietsch et alj l2005al and references therein). 
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Therefore, we identify the supersoft X-ray transient in the glob- 
ular cluster Bol 1 1 1 with the nova M31N 2007-06b reported by 
Shafter & Quimby (2007). This optical nova was first detected 
on 2007 June 19.38 and was detected in X-rays the first time 
141 days later on 2007 November 07.64. SSI was not visible on 
2007 July 22.55, which is 33 days after the first detection in the 
optical (see TableHJ. 



2.2. Supersoft source in Bol 194 

A second supersoft source (SS2) was detected in the same 
Chandra HRC-I observ ation as SSI (ObsI D 8526) as a new, 
bright X-ray transient dHaberl et al.l l2007ab . In the following 
observation (ObsID 8527) SS2 was again detected with sim- 
ilar brightness close to the edge of the HRC-I field of view 
(see Fig.|2]i. For the remaining three Chandra observations of 
the monitoring campaign the position of SS2 was outside the 
field of view due to the changing roll angle of the observations. 
Therefore, we used our Swift XRT ToO follow up observations 
(starting with ObsID 00031027001) to constrain the spectrum 
of SS2 and follow the X-ray light curve. The X-ray luminosi- 
ties given in Table|2] show that this source faded much faster 
than SSI. Therefore, no XMM-Newton observations of SS2 are 
available and we could use our XMM-Newton data only to com- 
pute upper limits. Since SS2 also has a large off-axis position in 
the HRC-I field, we used the Swift observation 00031027001 to 
compute the position of the source. We found the following coor- 
dinates: RA(J2000) 00:43:45.3, Dec(J2000) +41:06:08.15. With 
a lcr position error of 4 76 and a distance of 1 71 to the GC Bol 
194 dGalleti et alj|2004 00:42:45.20, +41:06:08.3) these coor- 
dinates are in good agreement with the position of Bol 194. See 
Fig.[3]for a comparison of optical (DSS POSS-II Red) and X-ray 
positions. 

The Swift XRT spectrum of SS2 (ObsID 00031027001) 
only shows photons with energies below 750 eV (see Fig.[7j). 
Therefore, we classified this source as a SSS. Since we have not 
many spectral counts for this source, we use Cash statistics for 
spectral modelling. The best spectral fit is an absorbed black- 
body with best fit values of kT = 74+g e V and N H = (l-0^) 
xlO 21 cm 2 . The best fit parameters differ fro m the preliminary 
data analysis published in lHaberl et alJd2007b|). where we used a 
differ ent abundance table in XSPEC (angr; Anders & Grevessel 



119891) and x 1 statistics. Using the best fit values we computed 
the X-ray luminosity of SS2 to be (6.4 + 1.0) xlO 37 erg s~\ the 
bolometric luminosity to be 10^ xlO 37 erg s _1 and therefore the 
radius of the WD to be 5*1 1 xlO 8 cm. Unfortunately, the low 
count Swift XRT spectrum does not provide enough degrees of 
freedom to test the WD atmosphere model that we fitted to the 
XMM-Newton EPIC P N spectra of SS I. 

After our discovery. iGalachee t al. (2007) reported that they 
found SS2 to be present in a 5ks Chandra ACIS-S observation 
taken on 2007-1 1-03. 18 UT, which is four days earlier than our 
first Chandra HRC-I detection. They confirm our classification 
of SS2 as a SSS and further report that the source is not visi- 
ble anymore in a 5ks Chandra ACIS-S observation on 2007-11- 
27.16 UT (~ 24 days after the previous ACIS observation) with 
an 95 % upper detection limit of ~ 6 xl0~ 4 ct s -1 . Note, that 
this observation was taken just three days after our Swift ob- 
servatio n on 2007-1 1-24.34 UT Ad ditional Swift XRT observa- 
tions bv lKong & Di Stefanol d2007al) showed a re-brightening of 
the source on 2007-12-02.64. However, the source faded quickly 
again and was not found in follow-up Swift observations and our 
XMM-Newton monitoring data. We re-analysed the Swift and 



% o.oi 



U 5x10-' 



0.2 



0.5 

Channel Energy (keV) 



Fig. 7. Swift XRT spectrum of SS2 in Bol 194 from observa- 
tion 00031027001 (7.3 ks exposure time) fitted with an absorbed 
blackbody fit (solid line). 



Chandra HRC-I data, computed upper limits for the XMM- 
Newton observations and present all available data on SS2 in 
Table[2] We took the count rates fo r detections and upper lim- 
its from lKong & Di Stefanol d2007al) to compute luminosities in 
XSPEC using our spectral model. For the Swift observations 
00031027004/5 we give 95% upper limits on the SS2 luminos- 
ity. These values are comparable to actual luminosities obtained 
for earlier Swift observations, an effect which is due to longer 
effective source exposure times for the earlier observations. X- 
ray luminosities for the three instruments were computed using 
XSPEC and the blackbody fit inferred from the Swift spectrum. 

Analogous to SSI and its possible connection to 
M31N 2007-06b in Bol 111, the supersoft spectrum and 
the transient light curve could indicate that SS2 is the X-ray 
counterpart of a recent optical nova in the GC Bol 194. The 
following section describes the optical data analysis and the 
constraints that we can put on a possible counterpart nova in 
Bol 194. 



3. Search for an optical nova counterpart of the 
supersoft source in Bol 194 

As there was no optical nova reported in Bol 194 we searched the 
recent optical data available to us for indications of a nova out- 
burst in the GC. Note, that the constant background light from 
the GC makes it difficult to detect nova outbursts. However, dif- 
ference imaging is a successful method in this context and was 
used to analyse a big part of our optical data set. In the follow- 
ing we describe the observational setup and the analysis of the 
optical data. 

3.1. Optical observations 

The majority of our optical data wa s obtained in t he context 
of the Texas Supernova Search (TSS; IOuimbvll2006l) . The TSS 
employed the 0.45-m ROTSE-IIIb telescope at the McDonald 
observatory in Texas, and the data presented here are supple- 
mented by its twin, the ROTSE-IIId telescope, which is lo- 
cated at the Turkish National Observatory at Bakirlitepe, Turkey . 
The ROTSE-III system is described in lAkerlof et all d2003l) . 
The 1°85 x 1°85 field of view covered by the unfiltered, 2kx2k 
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Marconi CCD encompasses practically all of M 3 1 's light in a 
single exposure (see Fig. and additional overlapping fields 
to the northeast and southwest add to the haul. Beginning in 
November 2004, these 3 fields were imaged several times nightly 
as weather and season allowed. 

Our second monitoring program for optical novae uses 
the robotic 60 cm telescope with an E2V CCD (2kx2k) 
Livermore Optical Transient Imaging System (Super-LOTIS, 
IWilliams et al.ll2008l) located at Steward Observatory, Kitt Peak, 
Arizona, USA. Starting in October 2007 the telescope was used 
in every good night to monitor the bulge of M 3 1 . Using four 
Super-LOTIS fields (field of view: 17' x 17') these observations 
cover an area of ~ 34' x 34' centered on the core of M 3 1 with 
a pixel scale of 0.496"/pixel and a typical limiting magnitude of 
19.0 mag, using a Johnson R filter. The reduction of the data is 
done by a semi-automatic routine and the astrometric and photo- 
metric cahbratiojn^esjhe M 31 part of the Local Group Survey 
(LGS, iMassev et al.ll2006l) . Typical values of Super-LOTIS as- 
trometric and photometric lcr accuracies are 0'.'25 and 0.25 
mag, respectively, averaged over the whole magnitude range. 
The southeast Super-LOTIS field covers Bol 194 and the astro- 
metric and photometric lcr accuracies for this particular object 
(~ 16.8 mag) are O'.'l 1 and 0.06 mag, respectively. 

Additionally, in this research our optical data set is supple- 
mented by archival data from K. Hornoch obtained at telescopes 
in Lelekovice (Newtonian focus of 350/1658 mm telescope, 
CCD camera G2CCD-1600, Kron-Cousins R filter, 1.12"/pixel, 
FOV 28.7'xl9.1') and Ondfejov (primary focus of 650/2342 mm 
telescope, CCD camera G2CCD-3200, Kron-Cousins R filter, 
1.20"/pixel, FOV 21.8'xl4.7'). Standard reduction procedures 
for raw CCD images were applied (dark and bias subtraction 
and flat-field correction) using SIMS0and MunipaclQ programs. 
Reduced images of the same series were co-added to improve the 
S/N ratio (total exposure time varied from 600s up to 1 800s). The 
gradient of the galaxy background of co-added images was flat- 
tened by the spatial median filter using SIMS. These processed 
images were used for aperture photometry, carried out in GAIAB 
Relative photometry was done using brighter field stars which 
were calibrated using standard Landolt fields. The lcr measure- 
ment uncertainties were low (~ 0.03 mag), thanks to the long 
exposure times, the brightness of Bol 194, and its location far 
from the high surface brightness levels found near the center of 
M31. 

3.2. Optical data analysis 

We searched the optical data of Bol 194 obtained by the moni- 
toring programs described above for a significant optical excess 
that could indicate a nova outburst in the GC. The significance 
for the Super-LOTIS data and for Hornoch's data is defined by 
applying a 3<x clipping to the light curves of Bol 194 and judging 
every two neighboured data points that lie outside the final +3cr 
range as possible indications of a nova outburst. 

For the ROTSE-III data, we used the PSF-matched im- 
age subtr action code developed by the Supernova Cosmology 
Project (Perl mutter et al.l 1 19991) to search for residual light 
from an optical nova. With this method th e outburst of nova 
M31N 2007-06b in Bol 1 1 1 was detected dShafter & Ouimbvl 
|2007|) . We first constructed a deep reference image by co-adding 
100 ROTSE-IIIb images obtained between December 2004 and 

s http://ccd.mii.cz/ 

6 http://munipack.astronomy.cz/ 

7 http : //www . starlink . rl . ac . uk/gaia 



June 2006. We then convolved this image to match the PSF of 
each ROTSE-III image, subtracted off this template light, and 
searched for any point sources coincident with Bol 194. We 
attempted to measure any resid uals at this lo cation with the 
D AOPHOT PSF-fit ting routines dStetsonlll987t ported to IDL 
bv lLandsmanll 19891) . 

These procedures did not lead to the discovery of an opti- 
cal nova counterpart of SS2 in Bol 194. However, our optical 
monitoring data of the M 3 1 bulge region allows us to put strong 
constraints on the outburst date of any nova that may have oc- 
curred in Bol 194. Figure [8] includes all our optical data of Bol 
194 from November 2004 until 2007-1 1-08, which is 5 days af- 
ter the first detection of SS2 in X-rays on 2007-1 1-03. We show 
the minimum detectable magnitude of a possible optical nova 
occurring in Bol 194 during this period. 

To calculate the limiting magnitudes of the ROTSE-III obser- 
vations, we measured the noise on the subtracted frames in an- 
nuli centered on the location of Bol 194. The limits reported here 
correspond to the flux required of a point source to be detected 
at the Act level. The magnitude scale was calibra ted against the 
USNO-B1.0 R2 measurements dMonet et al.ll2003l) . 

In contrast to this procedure, the limiting magnitude for the 
Super-LOTIS data and for Hornoch's data is computed as fol- 
lows. At the distance of M 31 the light of a nova in a GC would 
blend with the light of the GC itself. Therefore, we have to take 
into account the intrinsic magnitude of the cluster (~ 16.8 mag 
in the Super-LOTIS data) and compute the resulting magnitude 
using 

R, = -2.5/ O g 10 (10- a4R - + 10- a4s »), 

with R,, R gc , and R n being the total magnitude, the magnitude 
of the GC, and the magnitude of the nova, respectively. Thus, 
for the Super-LOTIS data and for Hornoch's data, the minimum 
detectable magnitude R/ of a nova is defined as the magnitude 
that leads to a significantly brighter/?,. The significance criterion 
used here is: 

R, < R gc - 3cr Rs( , 

with 3cr Sg being the 3<x photometric standard error for Bol 194 
(~ 0.16 mag for the Super-LOTIS data reduction and ~ 0.09 
mag for Hornoch's data). These measurement errors are mean 
values derived from our light curves of Bol 194. Since we used 
a 3<x clipping method to search for a nova outburst in Bol 194, 
we assume that the mean photometric errors given for both in- 
struments would correspond to the final 3cr range in the case of 
a nova outburst-modified light curve. Therefore, we compute R/ 
as follows: 

Ri = -2.5 Zog 10 (10-°' 4(R ^ 3<r v> _ 10-°- 4s *0- 

In order to judge the time coverage of our data set, we sim- 
ulated the outburst of novae with different peak magnitudes and 
decay times in Bol 194 on any day between 2004 October and 
2007 November and checked whether we would have detected 
the outburst according to the limiting magnitudes, or not. We 
simulated novae with peak magnitudes of 15.5, 16.0, 16.9, 18.0, 
18.5 and 19.0 mag in the R band and decay times of 6, 15, 
29, 63, 102 and 164 days, respectively. The peak magnitude of 
nova M 3 IN 2007-06b in Bol 1 1 1 as given in lShafter & Ouimbvl 
d2007l) is 16.9 mag. The decay time ti is defined as the time in 
days the nova luminosity needs to drop 2 mag below peak lu- 
minosity ( Pavne-Ga poschkin 1964). To compute this value, we 
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used the maximum magnitu de versus rate of de cline (MMRD) 
relationship given bv lDella Valle & Livid (fl995h . for t 2 < 50 d: 

M v ,max = -7.92 - 0.81 arctan J^__lpflJ > 

and bv lCohenl (fl988h . for f 2 > 50 d: 



Mv 



= 2.41 logf 2 - 10.70. 



We us ed this c ombination of the two different MMRDs ac- 
cording to IShafterl d 19971) i n order to compensate for t he break- 
down of the MMRD from iDella Valle & Livid d 1995b at large 
values of f 2 . Both MMRDs give the relation of the decay time of 
a nova to its maximum magnitude in the V band. The color re- 
lation for novae at maximum lig ht is (B - V)q = 0.23 + 0.06 
( Ivan den Bergh & Younger! [l98 7). This fact implies for novae 
at maximum light a spec tral type clo se to F0 and therefore a 
(V - R) ~ 0.30 (see e.g. lAllenlfl 9761) which we used to match 
both MMRDs to our R band data. 

For our simulations we u sed a M 3 1 extinction free distance 
modulus of /io = 24.38 mag ([Freedman et al. 2001). We applied 
a reddening towards M 31 of E(B - V) = 0.062 dSchlegel et alj 
1998) and a foreground extinction in the R band of A^ = 0.17 
mag, obtained via the NASA Extragalactic Database (NED). We 
estimated the internal absorption of M 31 to be A' R ~ 0.1 1 mag, 
using A' ~ 0.20 mag (Capa ccioli et al.ll 989b and assuming that 
A pg ~ A B . 

The results of these simulations are shown in the top panel of 
Fig.[8j with the upper parts of the curves indicating the periods 
over which each simulated nova would have been detected. The 
figure visualises the fact that our monitoring is efficient in the 
sense that the only periods in which we would not have been able 
to detect the most novae, are the periods when M 3 1 is not ob- 
servable due to its annual visibility. These periods, computed for 
a nova with parameters of M31N 2007-06b, and the estimated 
constraints on the delay time between optical outburst and detec- 
tion in X-rays are shown in Table[4] Note also, that a bright nova 
(peak luminosity = 15.5 mag) would have been missed more 
likely than a fainter nova. This is due to the faster decline of lu- 
minosity for the brighter novae and good detection limits of our 
observations down to ~ 19 mag. 



Table 4. Optical constraints on a nova in Bol 194 



Datel [UT]° 


Date2 [UT]" 


RJD1 [d]" 


RJD2 [df 


Delay [d] f 


2005-02-15 


2005-05-12 


3417 


3503 


990 - 904 


2006-02-19 


2006-05-09 


3786 


3865 


621 - 542 


2007-02-18 


2007-05-20 


4150 


4241 


257 - 166 



Notes: ": Start date (datel) and end date (date2) of the non- 
detection periods 
b : The same as in a but in RJD = JD - 2 450 000 
c : Periods of possible delay time between nova outburst 
and first detection in X-rays 



4. Discussion 

4.1. Summary of observational data 

We present in this paper the first two SSSs found in M 31 GCs. 
The source SSI in Bol 111 can be identified with M31N 2007- 
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Fig. 8. Top: Limiting magnitudes of all optical data since 2004 
November. Bottom: simulated detection (upper points) or non 
detection (lower points) for novae with indicated peak magni- 
tudes. In both panels the vertical dashed line shows the date of 
the first detection of SS2 in X-rays. 



06b, the very first nova found in a M 31 GC in June 2007 
dShafter & Ouimbvll2007l) . After 137 days (November 2007) we 
detected supersoft X-ray emission from this source and it was 
still active in our last XMM-Newton observation in July 2008. 
The supersoft spectrum and the transient nature of SS2 in Bol 
194 could also indicate a recent nova outburst. However, we 
could not find evidence for an nova counterpart in optical mon- 
itoring data of this GC. We give constraints on a possible nova 
outburst in Bol 194, based on the detection limits of our opti- 
cal observations and simulated nova outbursts. The simulations 
show that most novae would have been detected in Bol 194, ex- 
cept in the times when M 3 1 could not be observed due to it's an- 
nual visibility. We computed the time spans in which a nova like 
M31N 2007-06b in Bol 194 could have had an unobserved out- 
burst. The associated delay times, given in Tabled between op- 
tical outburst and X-ray switch-on are physically feasible, since 
similar dela ys have been observe d in other supersoft nova coun- 
terparts (see lPietsch et al.il2007bj) . Therefore, a nova outburst in 
Bol 194, that is connected to SS2, cannot be ruled out. 

4.2. Globular cluster nova rate - the optical point of view 

The statistics of nova e in GCs of other galaxies is still poor. 
Tomaney et alj (1 19921) conducted a search in Ha of over 200 
M 31 GCs for nova eruptions and found nothing over an ef- 
fective survey time of one year. This results in an upper limit 
for the M 31 GC nova rate of 0.005 novae yr~' GC" 1 . Anoth er 



search f or novae in 54 M 31 GCs was done by Ciardullo et al. 
(1990b) based on the M 31 Ha survey data of ICiardullo et al. 
(119871 11990a). Over a mean effective survey time of approxi- 
mately two years no indications for a nova outburst in one of the 
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GCs was found. From one nova fo und in a GC of the giant ellip- 
tical galaxy M 87 lShara etail d2004 derive a rate of 0.004 novae 
yr 1 GCT 1 for the whole system of 1057 known M 87 GCs (4.2 
novae yr -1 ). 

There are about 500 GCs known that belong to M 31 (482 
confirmed GCs, ac cording to the Revised Bologna Cata loguj^ 
(V.3.5, March 2008: lGalleti et al.l2004ll2005ll2006ll2007M This 
corresponds approximately to a stellar mass of about 5 xlO 9 
M©. Given a rate of 2.2 novae yr -1 (10 10 Mq) _i , which is 
observed in old systems like elli p tical g alaxies (data derived 
from Table 3 of iDella Valle et al.1 (11994 and Equation (1) of 
Man nucci et al.l (120051) ). we derive a rate of about 1 . 1 novae yr -1 
for the M 31 GC system, which corresponds to about 0.002 no- 
vae yr -1 GC -1 . This figure is comparable (within a factor of 2) 
wi th the observed ra te (0.004) of the M 87 GC system obtained 
by IShara et al. I (12004 and sugge sts that the discovery of only 
one nova in the M 3 1 GC system (IShafter & Ouimbvl20 07) over 
about a century of obse rvations may be the re sult of a strong 
observational bias (see [Shafter & O uimbv 20071) . 

4.3. Globular cluster nova rate - the X-ray point of view 

While the detection of CNe in GC in the optical is strongly ham- 
pered by the light from the GC itself, the detection of supersoft 
emission from a hydrogen-burning post-nova atmosphere is not 
affected by this. However, the presence of other X-ray sources 
in a GC, such as LMXBs, may create a similar problem for the 
detection of a SSS in that GC. In the following, we discuss the 
impact of our results on the rate of novae for the M 31 GCs. 

If we assume that both SSI and SS2 are post-novae and 
assume the duration of the SSS phase to be of the order of 
1 yr then we find a nova rate of 0.015 novae yr -1 GC" 1 for 
the ~ 130 GCs from iGalleti et all d2007l) in the field of view 
of our Chandra HRC-I observations (see also Fig. [2] for the 
Chandra field). The duration of the SSS phase of novae is 
a critical parameter for a correct estimate of the nova rate. 
Although there are some n ovae known to have m uch longer SSS 
phases than 1 yr (se e e.g. lOgelman et alj[l993h it is likely (see 
Piets ch et alj|20 07b) that novae with short SSS phases were just 
not found in the past due to selection effects. Also note, that the 
SSS phase of SS2 lasted only < 4 months and therefore the as- 
sumed duration of 1 yr may still be a quite conservative estimate. 

The nova rate computed using the X-ray data on the SSSs 
is about a factor 10 larger (after assuming the duration of the 
SSS phase to be 1 yr) than the nova rate estimated by assum- 
ing that GCs produce novae with a rate comparable to the old 
stellar populations in Ellipticals. The match between the rates 
in X-ray and in optical is obtained for a duration of the SSS 
phase o f about a decade. Su ch a long time has been observed 
indeed dO gelman et al. 1993), but it does not seem to be the rule 
dPietsch et al.ll2007bl) . 

The results based on the small size of our sample of GC SSSs 
are clearly dominated by Poissonian statistics. In this case, the 
discovery of two GC novae would still be consistent with a nova 
rate of ~ 0.002 novae yr -1 GC 1 (95% confidence level) and ~ 
0.004 novae yr -1 GC 1 (99.87% confidence level), respectively. 
Both values are below the upper li mit of 0.005 novae yr~' GC -1 
estimated bv lTomanev et al.l (Il992l) from optical observations. 

Furthermore, we want to emphasize the fact that not all op- 
tical novae in the Galaxy and in M 31 were detected as SSSs. 
Novae with supersoft X-ray emission are a subset, with yet un- 
known size, of the whole nova population. The nova rates de- 

8 http://www.bo.astro.it/M31/ 



duced from X-ray observations are therefore lower limits on the 
actual nova rates. 

4.4. Conclusions 

The nova rate in the M 3 1 GC system may be larger by about one 
order of magnitude than expected in old stellar systems like gi- 
ant ellipticals. This fact implies the existence of some other (dy- 
namical) mechanism which acts inside the GCs and is capable 
to raise the fraction of CN systems naturally produced by stellar 
evolution. It is well known that LMXBs are highly ove rabundant 
in GC s, with respect to the rest of a galaxy (see e.g. iFan et all 
120051 and references therein). This is explained by dynamical ef- 
fects like tidal captures of low- mass main-seq uence stars by neu- 
tron stars (see also lOarkll 19751: iFabian et aT]|197 5). One would 
expect similar effects for WDs dHut & Verbuntlll983l) . but al- 
though large numbers of cataclysmic variables (CVs) have been 
found in Galac tic GCs, the num ber of novae detected in GCs is 
still small (see Sharaetal . 2004, and references therein). 

Further, we note that novae with short SSS states seem to 
be an importa nt contributor to t he SSS population of galax- 
ies (see also lPietsch et al.ll2007bl). XMM - Newton observations 
of M 31 in the past dStiele et al l 120081; iPietsch et alj l2005bt 
iTrudolvubov & Priedhorsky 2004) did not detect SSS with po- 
sitions that agree with GC positions. This might in part be due 
to the fact that the individual observations of these monitorings 
were separated by half a year or more. The discoveries of SSI 
and SS2 are the first results of a new monitoring strategy, that 
uses XMM-Newton and Chandra observations that are separated 
just by 10 days. Future observations that follow this strategy may 
increase the statistics for SSS in M 31 GCs. If future surveys of 
M 31 will be able to find between six to eight new GC novae 
per year, the existence of a nova rate excess, of the size reported 
above, will be proven at 95% and 99.87% confidence level, re- 
spectively. 

Finally, given the possible li nk between SSS and type 

la supernovae (SNe-I a) (e.g. IDella Valle & Livid 119941: 

iDi Stefano & Rappaportlll994 if a high rate for SSS in GCs 
will be confirmed by future surveys, it may not be inconceivable 
to expect detection of SNe-Ia in GC systems around giant 
ellipticals or bulge dominated galaxies. 
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